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Abstract. - We report a numerical analysis of the Anderson transition in a quantum-chaotic 
system, the quasiperiodic kicked rotor with three incommensurate frequencies. It is shown that this 
dynamical system exhibits the same critical phenomena as the truly random 3D- Anderson model. 
By taking proper account of systematic corrections to one-parameter scaling, the universality of 
the critical exponent is demonstrated. Our result v — 1.59 ±0.01 is in perfect agreement with the 
value found for the Anderson model. 
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Introduction. — It is now widely acknowledged that 
the classical diffusive behavior of non-interacting electrons 
in a disordered potential can be stopped by non-trivial 
interference effects [1]. This puzzling phenomenon, An- 
derson localization, constitutes one strong evidence of the 
very difference between quantum and classical dynamics 
of complex systems. In the field of quantum chaos, a sim- 
ilar phenomenon is acknowledged, namely dynamical lo- 
calization. In the kicked rotor, a paradigmatic system of 
quantum chaos, quantum mechanical interference tend to 
suppress the classical chaotic diffusive dynamics. The dis- 
covery of the parallel between dynamical localization and 
Anderson localization originated from the mapping of the 
kicked rotor to the quasirandom ID Anderson model [2], 
In Ref. [3] it was demonstrated that the kicked rotor in the 
dynamical localization regime could be modeled by ran- 
dom band matrices; the latter have been reduced to a ID 
nonlinear a model [4] similar to those employed in the lo- 
calization theory [5] . In Ref. [6] the direct correspondence 
between the kicked rotor and the diffusive supersymmet- 
ric nonlinear a model was demonstrated. In the localized 
regime, the kicked rotor exactly mimics the behavior of 
disordered electronic conductors. 

There has been much efforts to observe Anderson local- 
ization in 3D experimentally. However, due to stray ef- 
fects like interaction, decoherence or absorption, very few 
attempts have been successful [7]. In a slightly different 



context, Anderson localization of acoustic [8] and electro- 
magnetic [9-12] waves has been experimentally observed. 
The experimental realization of the kicked rotor with laser- 
cooled atoms interacting with a pulsed standing wave al- 
lowed for the first experimental observation of Anderson 
localization in ID with atomic matter waves [13]. One step 
further is to observe the well-known Anderson transition 
with this type of system, i.e. the disorder induced metal- 
insulator transition predicted for non-interacting electrons 
subjected to a 3D disordered potential. Different general- 
izations of the kicked rotor have been theoretically consid- 
ered as analogs of the 3D- Anderson model [14]. Here, we 
focus on the convenient three-incommensurate- frequencies 
generalization introduced in Ref. [15]. Very recently an ex- 
periment based on this system has fully characterized the 
Anderson metal-insulator transition [16]. Indeed, a care- 
ful analysis of the scaling properties of the dynamics has 
resulted in the first experimental determination of the lo- 
calization length critical exponent v of the Anderson tran- 
sition. The value found v = 1.4 ± 0.3 is compatible with 
the precedent numerical determination of v = 1.57 ± 0.02 
for the true-random 3D- Anderson model [17]. 

At this stage, the equivalence between the quasiperi- 
odic kicked rotor [16] and 3D-disordered conductors still 
has the status of a conjecture (see [18]). A rigorous an- 
swer to the question whether this dynamical system ex- 
hibits the same critical phenomena - i.e. belongs to the 
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same university class - as the true 3D-Anderson model 
has not yet been given. Can a simple three-frequency dy- 
namical system exactly mimic the critical behavior of 3D 
disordered electronic conductors? In this Letter, we show 
that the answer is positive. This is done by carrying out 
a very precise numerical study of the critical behavior of 
the quasiperiodic kicked rotor with the same care as in the 
most sophisticated investigations of the critical behavior 
of the true 3D- Anderson model [17,19,20]. The fact that 
both models give the same localization length critical ex- 
ponent v within comparable (small) uncertainties implies 
that they belong to the same universality class (orthogo- 
nal) [21]. 

The quasiperiodic kicked rotor. — The quasiperi- 
odic kicked rotor we consider is a three-incommensurate- 
frequencies generalization of the kicked rotor: 



H qp = ^-+ K{t) cos eJ2 S (t- ") 



(1) 



obtained simply by modulating the amplitude of the 
standing wave pulses with a set of two new incommen- 
surate frequencies u 2 and u 3 of modulation: 

K.(t) = K [1 + ecos (uj 2 t + ip 2 ) cos (u 3 t + ip 3 )] . (2) 

Here we consider the case with an effective Planck constant 
h = —i[9,p\. For a standard rotor, is an angle defined 
modulo 2tt and the wavefunction thus has to be 2ir peri- 
odic. In the atomic realization of the kicked rotor [16], the 
Hamiltonian is still given by Eq. (1), with 8 extended in 
the (—00, +00) range. Using the Bloch theorem, one can 
restrict to 2-7T periodic functions, at the (cheap) price of 
including a constant quasi-momcntum. 

The dynamics of this quasiperiodic kicked rotor is iden- 
tical to the time-evolution of a 3D-kickcd rotor: 



Pi 

H 3 = — + uj 2 p 2 + u 3 p 3 

+K cos 6>i [1 + e cos 2 cos 3 ] S ( f ~ n ) ' ( 3 ) 



with an initial condition: 



1-3(01,02,03,* = 0) 

= * w (0i,t = O)5(0 2 - V2 )5(6 3 -<p 3 ) 



(4) 



where ^ qp {9, t = 0) is an arbitrary initial condition for the 
quasiperiodic kicked rotor. Note that dynamical localiza- 
tion takes place in momentum, not configuration, space. 
The initial state being perfectly localized in 02 and 03, it is 
entirely delocalized in the conjugate momenta p 2 and p 3 so 
that we will study transport along the pi direction, which 
is tantamount to measure the momentum distribution of 
the quasiperiodic kicked rotor \ 1 $ qp (p, t)\ 2 . The unusual 
linear dependence of H 3 with p 2 and p 3 does not prevent, 
for nonzero e the dynamics to be similarly diffusive or lo- 
calized along the 3 coordinates. 



The Hamiltonian H 3 is invariant under the transforma- 
tion T : t — > —t, 6 — > 6, p — > p, i.e. the time-reversal in 
the momentum representation, which is the relevant one 
for the dynamical localization (see [22,23]). In particular, 
the choice of a non-zero quasi-momentum and non-zero 
phases ip 2 or ip 3 docs not break time-reversal symmetry. 
The evolution of the states according to the Hamiltonian 
(1) is governed by the operator U [see below Eq. (5)], be- 
longing to the Circular Orthogonal Ensemble class, with 
the additional constraint at t = (4), this shows that the 
dynamical properties of the present quasiperiodic kicked 
rotor also belong to the orthogonal ensemble. 

It should also be noted that the 3D aspect comes from 
the fact that 3 frequencies are present in our dynami- 
cal system: the usual "momentum frequency" £ which 
is present in the standard kicked rotor (with e = 0), and 
the two additional time-frequencies lo 2 and lu 3 . By in- 
creasing the number of incommensurate frequencies, one 
should be able to tune the effective dimensionality of the 
system. This holds the promise of extending the study of 
the Anderson transition to higher dimensions. 

From a stroboscopic point of view, the quantum dy- 
namics of the 3D-kicked rotor Eq. (3) is determined by its 
evolution operator over one period: 



jj g— iK cos 0i cos 62 cos 63)/% 

Xe - i (Pl/ 2 +"2P2+^3P3)/S; 



(5) 



whose eigenstates form a basis set allowing to calculate 
the temporal evolution. These Floquet states \4>) are fully 
characterized by their quasienergy cu, defined modulo 2ir: 



U\<j> u ) = e-*"^) 



(6) 



Equivalence with a 3D- Anderson tight-binding model can 
be obtained by reformulating Eq. (6) for the Floquet states 
[2]: 

fm$m + ^m-r = -W <P m , (7) 

where m = (mi, m.2, ^3) and r label sites on a 3D lattice, 
and the <& m are simply related to the Fourier components 
of the Floquet state |^)- The on-site energy e m reads: 



tan ■{ - 
2 



mi 



+ aj 2 m 2 + u 3 m 3 



)]} 



(8) 



and the hopping amplitudes W r are coefficients 
of a threefold Fourier expansion of W(6) = 
tan [K cos 0i (1 + e cos 02 cos 03) / 2k] . 

When (U,oj 2 ,uj 3 ) is an incommensurate triplet of fre- 
quencies, each component of the triplet being in addition 
incommensurate with ir, the classical dynamics can be- 
come chaotic (for sufficiently large stochasticity parameter 
K > 2) with diffusive spreading in all m directions [15]. 
The pseudo-random character of the potential e m then 
gives the disorder in the model (7). 

We therefore expect to observe localization effects as 
predicted for the standard 3D- Anderson model. Localized 
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states would be observed if the disorder strength is large 
as compared to the hopping amplitudes. In the case of 
model (7), while the amplitude of disorder is fixed, the 
hopping amplitudes increase with K. Therefore, for large 
K a diffusive regime should be observed, while for small 
K (though sufficiently large as to prevent classical local- 
ization effects) a localized regime is expected. This was 
validated both numerically [15] and experimentally [16]. 

Finite-time scaling. — In the case of the kicked ro- 
tor, the initial classical chaotic diffusion is stopped only 
after a certain characteristic time, the localization time Tg 
which turns out to be roughly proportional to the localiza- 
tion "length" (characterizing the exponential localization 
in momentum space) [24]: T£ oc £. In 3D, the localiza- 
tion time scales as T£ ~ £ 3 . Thus for very large localiza- 
tion length, Te may exceed the largest time accessible (the 
maximum duration of a cold-atoms experiment is typi- 
cally 150 kicks [16] whereas numerical investigations can 
go up to 10 6 kicks). Consequently, it is vain to investigate 
the Anderson transition only from static properties of the 
quasiperiodic kicked rotor, such as the divergence of the 
localization length at criticality, which could be obtained 
only for t > Tg. 

Actually, these finite-time effects arc similar to finite- 
size effects for the study of classical or quantum phase 
transitions. We can generalize the usual scaling laws to 
cover our time-dependent problem. The single parameter 
scaling theory, successfully used for the standard (static) 
3D Anderson model [25-27], can be applied to analyze the 
dynamics (see [16]) and especially to determine the criti- 
cal properties of the Anderson transition, i.e. the critical 
exponents. 

The dynamics of the quasiperiodic kicked rotor is con- 
veniently studied by considering the time-evolution of the 
variance of the momentum distribution {p 2 ). We can make 
the following scaling hypothesis for this quantity [28]: 



(p 2 )=t k iF[(K-K c ), 



(9) 



with F a function characteristic of the transition which is 
to be determined and fci and &2 two exponents which can 
be constrained as explained in the following. K c is the 
critical value of the stochasticity parameter. Note that 
(p 2 ) is here averaged over several initial conditions (this 
corresponds to an average over the disorder). We must 
recover as t — > oo either a diffusive behavior (p 2 ) ~ Dt 
when K > K c or a localized dynamics (p 2 ) ~ £ 2 when K < 
K c . In the vicinity of the transition, the diffusion constant 
D ~ (K — K c ) s vanishes with the critical exponent s and 
the localization length £ ~ (K c — K)~ v diverges with the 
critical exponent v. Together with Wegner's law [29] s = v 
in three dimensions, this leads to k\ = 2/3 and k-i = 
1/32/ [30]. 

Therefore, to analyze the scaling of the dynamics, the 
behavior of the quantity A = {p 2 )t~ 2 / 3 should be investi- 
gated as a function of time and for various stochasticity 




Fig. 1: (Color online) Dynamics of the quasiperiodic kicked 
rotor in the vicinity of the critical regime. The rescaled quan- 
tity \n\(K,t) is plotted as a function of K for various values 
of time t ranging from t — 30 to t = 40000. The crossing of 
the different curves at a common point (K c ~ 6.4, In A c ~ 1.6) 
indicates the occurrence of the metal-insulator transition. The 
parameters are that of the set A (see Table 1). 



parameters. Such a study was undertaken in [16] and al- 
lowed for a successful experimental characterization of the 
Anderson transition. In order to tackle the problem of uni- 
versality of the critical behavior, we have to go one step 
further and study whether the critical exponent v changes 
when parameters such as £,W2,o>3 are modified. 

The following discussion is rather intricate but cannot 
be avoided in a rigorous study. Indeed, to reliably dis- 
tinguish the different universality classes of the Anderson 
transition requires a very precise determination of the crit- 
ical exponent; for instance, the value v = 1.43 ± 0.04 for 
the unitary symmetry class is very close to the one for 
the orthogonal symmetry class [19]. For a very accurate 
determination of v, possible systematic deviations to one- 
parameter scaling must be taken into account. 

Systematic corrections to scaling. Let us con- 
sider the scaling function T = In F in the vicinity of the 
critical point: 



In A = T 



(K - K c ) t x ^ v 



(10) 



One simple feature of the scaling hypothesis Eq. (10) is 
that when In A is plotted against K, the curves for different 
times t should intersect at a common point (JsT c ,lnA c ); 
this crossing, which indicates the occurrence of the metal- 
insulator transition, is clearly visible in Fig. 1. 

In practice, the data do not exactly follow Eq. (10). 
There are small systematic deviations to scaling. Here, 
we consider several ways in which such corrections can 
arise (see below for a detailed discussion of each correc- 
tion): (i) the presence of an irrelevant scaling variable, 
(ii) non-linear dependence of the scaling variables in the 
stochasticity parameter K and (iii) resonances due to the 
periods being well approximated by a ratio of small in- 
tegers, (iii) is specific to our three-frequency dynamical 
system, but (i) and (ii) were shown to play an important 
role in the standard Anderson model [17]. Note that the 
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Fig. 2: (Color online) Small oscillating corrections to the scal- 
ing behavior in the critical regime. The parameters are the 
following: k = 2.89, u 2 = 2ttv%, lu 3 = 2ttv / 13, K = 7.8 and 
e = 0.3. Color curves correspond to various choices of the 
phases if2 and tpz [see Eq. (1)] whereas the black curve results 
from a statistical average over different phases. The amplitudes 
of the quasi-resonant oscillations decrease as time goes on. Av- 
eraging over the phases kills the rapidly oscillating structures, 
while keeping all of the other dynamical properties unchanged. 



most important correction to scaling is a time dependence 
of In A at K — K c either due to (i) or (iii). 

The well-known deviations (i) and (ii) can be taken into 
account by adding to (10) extra terms [17]. (i) The scal- 
ing function T depends not only on the relevant scaling 
variable Xr (i-e. a function of K — K c ), but also on an 
irrelevant scaling variable i/j: 



In A = T 



(xrt 



l/Zv 



(ii) 



Since tj) is an irrelevant scaling parameter, its effects 
should vanish as t goes to infinity, thus y must be pos- 
itive, (ii) Nonlinearity in the relevant scaling variable Xr 
can be described by an expansion in terms K — K c up to 
order mR. 

To define a fitting function J-f, we can then make a 
Taylor expansion of the scaling function T up to order ur 
in Xrt 1 ^ u and nj in ipt~ y : 



T f {K,t) = J2 J2x?t m/3 "i> n t~ ny F„ 



(12) 



We now discuss the qualitative nature of the correction 
due to the presence of an irrelevant scaling variable (i) 
as compared to the effect of resonances (iii). In case (i), 
lnA c = \nA(K = K c ) will shift in a monotonous way as 
time increases and converge to its thermodynamic limit 
lnA c (t = 00) [lnA c (i) = Ta.o + i _y J"o,i m the linear 
regime (to_r = tir = nj = 1)]. In the case of our three- 
frequency dynamical system, the data do not always fit 
such a monotonous evolution model for lnA c (i). Indeed, 
for a generic choice of incommensurate periods, the data 
are found to oscillate around their transient anomalous 
diffusive dynamics, see Fig. 2. 

We infer that such an oscillating correction to scaling 
arises in an resonant way: when the frequencies can be 
approximately related by a simple linear combination (i.e. 



involving small integers), a resonance occurs. This is a 
rather common phenomenon in multi-frequency dynamical 
systems. The phase and the amplitude of the oscillations 
in Fig. 2 depend on the choice of the phases ip2 and (^3 of 
the time-modulation [see Eq. (2)], i.e. on the initial state 
in Eq. (4). From the point of view of the Anderson-like 
model Eq. (7), resonances can be interpreted as correla- 
tions in the disordered potential. Hence, to perform the 
standard scaling analysis devised for the Anderson model 
with uncorrelated disorder [17], we shall retain data only 
for sufficiently long times and average over different ini- 
tial conditions, i.e. different quasi-momenta and phases 
ip 2 and </? 3 . 

We computed In A for times up to t — 10 6 kicks with 
an accuracy of 0.15%. To achieve this accuracy more than 
1000 initial conditions are required. By discarding data 
for times t < 1000, an analysis of the data that neglects 
oscillating corrections becomes possible. To analyze data 
over the full range of times t G [lO 3 , 10 6 ] , we fit the model 
Eq. (12) to the data. Note that the inclusion of the cor- 
rections (i) and (ii) in Eq. (12) leads to a rapid increase in 
the number of fitting parameters. That is why high qual- 
ity data with a wide range of variation of i 1 / 3 are needed 
if meaningful fits are to be obtained. 

The most likely fit is determined by minimizing the x 2 
statistics measuring the deviation, due to the numerical 
uncertainties, between the model and data. Some typi- 
cal numerical data and the associated fit are displayed in 
Fig. 3. To exhibit scaling, we subtract the corrections due 
to the irrelevant scaling variable [17] obtaining the cor- 
rected quantity lnA s . As seen in Fig. 4, all data collapse 
almost perfectly on the scaling function deduced from the 
model Tf. 

No significant deviation of the scaling function from the 
fit is observed. The goodness of fit Q has been determined 
using the x 2 distribution with Nd — N p degrees of freedom 
where Nd is the total number of data we used to fit the 
model [31] and N p is the total number of fitting parame- 
ters (see Table 1). The confidence intervals for the fitted 
parameters were estimated using the bootstrap method 
which yields Monte-Carlo estimates of the errors in the 
fitted parameters [31]. In the following, we shall present 
the fitted parameters along with the corresponding 68.2% 
confidence intervals (one standard deviation). 

Universality. — A key property of the Anderson tran- 
sition is that it is a continuous (i.e. second order) phase- 
transition in which the critical behavior can be described 
in a framework of universality classes [5] . This means that 
the critical exponents should not be sensitive to the mi- 
croscopic details of the disordered potential but should 
depend only on the underlying symmetries (e.g. time- 
reversal symmetry). We present here new material that 
allows us to numerically prove that this is indeed the case 
for the quasiperiodic kicked rotor, thus strengthening the 
fact that this system is in the same universality class than 
the truly random Anderson transition. 
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Fig. 3: (Color online) In A as a function of time t. The curves 
are fits of the data according to the model Tj [see Eq. (12)] 
with rifi — 3, niR = 2, m = 1. The parameters are that of the 
set T> (see Table 1). The inset shows the deviations of the data 
(corresponding to K = 7.9) to the most likely fit, showing no 
statistically significant deviation. 



We have carried on a detailed study of four cases char- 
acterized by different set of parameters, see Tabic 1): A 
and B are both optimal set of parameters for experimental 
studies (see [16]), while C and T> are rather for theoreti- 
cal/numerical considerations. C is a first step towards an 
ideal choice of parameters: 1^2/^3 — f] where 77 is the sil- 
ver number (see below), and the continuous fraction of 
H/tt is constituted of small integers (to prevent the sys- 
tem to be close to a resonance). V should be a "best 
choice" of parameters if we seek the least correlations in 
the disorder Eq. (8). ft is such that ft, u^, w 3 and ir are a 
"most incommensurate" quadruplet of numbers. As sug- 
gested in [14], we set k — a, 0J2 = a/rj and W3 = a/77 2 
where 77 = 1.324717... (the silver number which general- 
izes the golden number for triplet instead of pair of incom- 
mensurate numbers [32]) is the real root of the equation 
77 s — 77 — 1 = 0, and a = 3.5399... is such that the continuous 
fractions of k/ir — ct/ir, 102/^ = a/rjn and 0^3/71 = a/rj 2 ir 
are constituted of small (< 9) integers. 

The details of the simulations and the types of fit used 
to analyze those sets are listed in Table 2. The estimated 
critical parameters and their confidence intervals are given 
in Table 3. Some typical data and scaling function are 
drawn in Figs. 3 and 4. 

The most important point to be drawn from Table 3 is 
that the estimates of the exponent v for the four different 
sets are in almost perfect agreement with each other and 
with the estimate of v based on numerical studies of the 
truly random Anderson model v = 1.57 ± 0.02 of the or- 
thogonal symmetry class [17]. Note also that in the case 
of the quasiperiodic kicked rotor, the critical stochasticity 
K c and lnA c depend on: (i) the anisotropy governed by 
the parameter e and (ii) k, UJ2 and 0J3. The dependence (i) 
of the critical disorder and critical In A on anisotropy is a 




ln(£A 1/3 ) 

Fig. 4: (Color online) lnA s , the data in Fig. 3 after subtrac- 
tion of corrections due to the irrelevant scaling variable, plotted 

versus In (c/t 1/3 ) where £, ~ \xr\ " and the scaling function 
deduced from the model J-f, Eq. (12) (black curve). The pa- 
rameters are that of the set T> (see Table 1). The best fit 
estimates of the critical stochasticity and the critical exponent 
are in this case: K c = 8.09 ± 0.01, lnA c = 1.64 ± 0.03 and 
v = 1.59 ± 0.01. 



typical feature of the Anderson transition in anisotropic 
solids [20]. The dependence (ii) follows from the relation 
between the initial "classical" diffusion constant and the 
parameters k, 102 and [24]. 

The Anderson transition with the quasiperiodic kicked 
rotor is a robust feature: indeed, the "naive" choices of 
parameters, set A and £>, lead to very clean critical be- 
haviors, as clean as for the sophisticated choices of pa- 
rameters, sets C and T>. Our experience is that for certain 
mutually incommensurate triplets (£, u>2, W3) systematic 
deviations to scaling (such as resonances) can occur for 
intermediate times, but eventually vanish. 

Conclusion. — A numerical analysis of the critical 
behavior of the quasiperiodic kicked rotor has shown that 
this quantum chaotic system exhibits the same critical 
phenomena as the truly random Anderson model, i.e. 
both systems belong to the same (orthogonal) universal- 
ity class [21]. By taking proper account of corrections to 
the scaling property around criticality, the universality of 
the critical exponent v for the quasiperiodic kicked rotor 
was demonstrated. The critical exponent v was deter- 
mined with an accuracy better or comparable to the one 
achieved in previous numerical studies of the 3D- Anderson 
model [17, 19, 20]. Our result v = 1.59 ± 0.01 is in perfect 
agreement with the value found for the orthogonal sym- 
metry class [17]. It is clearly distinct from the value found 
for the unitary class v = 1.43 ± 0.04 [19], and from the 
predictions v = 1 of the self-consistent theory of localiza- 
tion [34] and v — 1.5 of a recent ad-hoc refinement of the 
self-consistent theory [35]. 
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Table 1: The four sets of parameters considered: k, u>2 and CJ3 
control the microscopic details of the disorder, while e drives 
the anisotropy of the hopping amplitudes. In C, 1^2/^3 — r\ 
where r) is the silver number (see below), and the continuous 
fraction of k/n is constituted of small integers. T> is such that 
k — a, u>2 = a/n and W3 = a/if where r\ = 1.324717... is the 
silver number, i.e. the real root of the equation r; 3 — 77 — 1 = 0, 
and a is such that the continuous fractions of k/n = a/71", 
(jJ2/tt = a/r/iY and u>s/tv — a/rfn are constituted of small (< 9) 
integers [32]. 
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